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Abstract 

In a previous work, a perturbative approach to a class of Fokker-Planck equations, 
which have constant diffusion coefficients and small time-dependent drift coefficients, 
was developed by exploiting the close connection between the Fokker-Planck equations 
and the Schrddinger equations. In this work, we further explore the possibility of 
similarity solutions of such a class of Fokker-Planck equations. These solutions possess 
definite scaling behaviors, and are obtained by means of the so-called similarity method. 
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1 Introduction 



One of the basic equations employed to describe fluctuating macroscopic variables is the 
Fokker-Planck equation (FPE) PQ. It has found applications not only in physics, but also 
in other areas such as astrophysics [21 [3J, chemistry (H [5], biology [HI E], finance [SJ, etc. In 
view of its wide applicability, various methods of finding exact and approximate solutions of 
the FPE's have been developed [D El [10J HI 021 1131 E] . Most of the methods, however, are 
concerned only with FPE's with time-independent diffusion and drift coefficients. Generally, 
it is not easy to find solutions of FPE's with time-dependent diffusion and drift coefficients. 

One of the methods of solving FPE with time-independent diffusion and drift coefficients 
is to transform the FPE into a time-independent Schrodinger equation, and then solve the 
eigenvalue problem of the latter [U [15]. The transformation to the Schrodinger equation 
of a FPE eliminates the first order spatial derivative in the FPE and creates a Hermitian 
spatial differential operator. In [16] it was shown that this transformation can also be 
applied to transform FPE's, which have constant diffusion coefficients and time- dependent 
drift coefficients, into time-dependent Schrodinger equations. Based on this connection, a 
perturbative approach was developed and used to solve this special class of the FPE [TB] . 

The diffusion equation is known to be invariant under the scale transformation x — > ex, 
t — > e 2 t for any scale e. Thus one can find solutions, called the similarity solutions, of the 
diffusion equation with definite scaling behaviors by means of the so-called similarity method 
[T7] . Such solutions were not considered in [16]. It is the purpose of this work to investigate 
the possibility of similarity solutions of the FPE within the perturbative framework of [16]. 

2 Perturbative approach to Fokker-Planck equations 

The FPE of the probability density W(x, t) in (1 + l)-dimension is [TJ 



0_ 

Of 



W(x,t) 



- _ d_ 

dx 



D^(x,t) + ^D^(x,t) W(x,t), 



(1) 
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where D^\x,t) and D^(x,t) are the drift and diffusion coefficient respectively. The drift 
coefficient represents the external force acting on the particles and is usually expressed in 
terms of a drift potential U(x,t) according to D^'(x,t) = —dll(x,t)/dx. 

In this paper we will focus on FPE's with constant diffusion coefficients D^ 2 \x, t) = D > 
0. In this case, one can solve the FPE by exploiting the connection between the FPE and 
the Schrodinger equation [H [15] . Setting 

Mx,t) = exp{^^}w{x,t), (2) 

one can transfrom the FPE into the Schrodinger equation 

dip _ d 2 ip /U" U' 2 U \ 



where the prime and dot denote the derivatives with respect to x and t, respectively. Each 
solution of the time-dependent Schrodinger equation gives the corresponding solution of 
the FPE via Eq. (J2J). It is, however, generally difficult to solve the Schrodinger equation 
exactly when the potential is time-dependent. Various approximation schemes may have to 
be employed. 

The perturbative approach presented in [16] to solve the FPE containing a small param- 
eter in the drift potential is summarized as follows. Suppose U(x,t) = J2^Lo ^ n ^n{ x i t), 
where |A| <C 1 is a small parameter, and let ipi^x^t) = exp{S(x,t, X)/D}. Then from Q, we 
get the equation of S(x, t, A), 

S = DS" + S' 2 + U , (4) 
V=^U»-\u* + \iJ. (5) 

Substituting the series form of S(x,t, A) = Yl^Lo A n S' n (x,t) into (j4j) and collecting terms of 
the same order of A, one arrives at a set of differential equations that determine the functions 
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s = dss + s>* + (£us - \k + \u ) , (6) 

51 = DS'i + 2S' S[ + (f - l -U' Q U[ + It)i) , (7) 

5 2 = DS» + 2S' S> + S» + (f US - \u' U' 2 - \u? + \u 2 ) , (8) 



^ n = DS" + S 'kS'n-k + terms in U of the order A n , n > . (9) 

fc=0 

With So, Si, 5*2,... solved, we will have an approximate solution ip(x,t), and hence of W(x, t). 

In this paper, we shall be interested in the drift potential U(x, t) of the form U(x, t) = 
XUi(x,t), with all other U n — for n = and n > 1. From the probability density is 

^^eM^ele^E^A^-Al) ? (10) 

We assume the initial profile of W(x,t) of the unperturbed case (the diffusion case) to be 
the delta-function, i.e., W(x,t) — >■ <5(x) at t = as A — > 0. Then the solution So(x,t) is [16] 



S (x,*) = -§M47r£>*)-^, (11) 

leading to the probability density 

W (x, t) = = - 7 =1= e"^ . (12) 

Eq. f ll2p is the well-known solution of the diffusion equation with the delta function as the 
initial profile. It is evident that under the scale transformation x — > x — ex, t — > i = e 2 t, 
W (x,t) scales as W (x,t) = eW (x,i). 



3 Scaling behavior 

Motivated by the scaling form of Wo(x,t) in fTl2|) . we would like to seek similarity solutions 
of the perturbative FPE. To this end, let us first obtain the required scaling behaviors of 
W(x,t) and U(x,t) so that Eq. ([T]) is invariant under the scale transformation 

x = e a x, i = e b t, (13) 
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where the scaling exponents a and b are arbitrary real parameters. Suppose W(x,t) and 
U(x,t) scale as W(x,i) = e 1 W(x,t) and U(x,i) = e d U(x,t) with real exponent 7 and d, 
resepctively. Then Eq. ([1]) becomes 



For Eq.f fT4|) to have the same form as Eq. ([I]), one must have e~ 1+b = £-"i- d+2a = e~"i +2a . 
This implies that d = 0, b = 2a, and 7 arbitrary. It means that when the scaling exponent of 
t is twice that of x, and U(x,t) is scale- invariant, then the FPE is scale invariant, and thus 
admits solutions W(x, t) with arbitrary scaling exponent, which is dictated by the initial 
profile W(x,0). This result puts a constraint on the scaling property of the drift potential 



Applying these general results to Eq. ffTOl) . we see that since U(x,t) = XUi(x,t) is scale- 
invariant, so are all S n (x,t) (n > 1). The scaling of W is therefore solely determined by the 
solution exp{So/D} of the non-perturbed FPE. With So(x,t) given in (TTTj) . W(x, t) scales 
as W(x, t) = e a W(x, t) (see Eq. ffT2]) ). This is in accord with the scaling of the initial profile: 
W(x, 0) = 5(x) = e a 5(x). Our next step is to obtain the scale-invariant solutions for the 
Sn's. This is attained by the similarity method, which we describe below. 

4 Similarity method 

The similarity method is a very useful method for solving a partial differential equation 
which possesses proper scaling behavior. One advantage of the similarity method is to 
reduce the partial differential equation to an ordinary differential equation through some 
new independent variables (called similarity variables), which are certain combinations of 
the old independent variables. In our case, the 2nd order FPE can thus be transformed into 
an ordinary differential equation which may be easier to solve. We will illustrate the method 
by discussing the solution of below. The method applies to equations for other S n as well. 




(14) 



U(x,t). 



5 



With Uo(x,t) = and S given in (ITTj) . the equation for Si is 



Si = DS'{ - + f C/f + |c>i . (15) 



Let x and t transform according to Eq. ([13]) . Recall from previous discussions that f/i and 
S n (n > 1) are scale-invariant, i.e., U(x,t) = Ui(x,t) and S n (x,t) = S n (x,t) {n > 1). In 
terms of these scaled variables, Eq.( lT5|) becomes 

8i dx 2 i dx 2 dx 2 2 dt ' 1 J 

We demand that Si satisfies Eq. fn5|) and this requires b = 2a. To determine the form of the 
similarity solution, we have to first determine a new variable z(x, t) that is invariant under 
the scale transformation Eq. ([T5|) . From b = 2a, we can choose the scale- invariant similarity 
variable to be 

' = 7t ■ (17) 

Since both Si(x,t) and Ui(x,t) are scaling invariant, they must be functions of z only. 
Hence it is reasonable to express Si(x,t) and Ui(x,t) in terms of similarity variable z as 

S 1 (x,t)=y(z); U 1 (x,t)=u(z) . (18) 

Eq. (|T5i) can then be cast into an ordinary differntial equation 

y"(z) - ^ y\z) + \ u"(z) - ^ u\z) = . (19) 

Eq. (1191) is the most general equation for S% with an arbitrary scale-invariant potential 
U (x, t) = XUi(x, t). Once y(z) is solved, the similarity solution Si(x, t) of Eq. (fT5|) is obtained 
by putting z — xj \fi. To illustrate the procedure, we shall solve the case of a class of simple 
scale-invariant potentials. 



5 Ui{x,t) = iix p t q 

For definiteness we consider the simplest form that may have the required scaling property, 
namely, Ui(x,t) = [ix p t q where /i, p and q are arbitrary real parameters. Requiring that 
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Ui be scale-invariant under the transformation (TT5|) with b = 2a, we must have q = —p/2. 
From ( lT8j) . we have 

u(z) = U x (x, t)=fix p t q = fiz p . (20) 

Eq.f lTgj) is then expressed as 

y" -—y' = -^p(p - 1) z p - 2 + ^z p . (21) 
y 2D y 2 FKF 'AD y ' 

The general solution of Eq. fT2T]) is 

y(z)=a J' exp + z*. (22) 

In order for the probability density W(x, t) to be normalizable, we require that y(z) be finite 
as z — > ±oo. As such, ao has to be set to zero, as the exponential term is divergent in the 
domain of z G (— oo, oo). Thus the acceptable solution is 

y(z) = «i - \z p . (23) 

Putting z = x/y/t back into (|23|) gives the solution of the first order perturbation equa- 
tion: 

S 1 (x,t) = a 1 -^(Jj=y . (24) 
With Si(x,t) given in (|2~4"|) . the equation for S 2 (x, t), i.e. Eq.([H]), is 

S 2 = DS'i - jS> 2 . (25) 
By the procedure described in Sect. 4, this equation can be reduced to 

v" -^y' = °- ( 26 ) 

As discussed before, the finite solution is y(z) = a 2 = constant. Hence 5* 2 is given by 
S 2 (x,t) = a 2 . 

With S\ and S2 given above, the equations for S n {n > 3) turn out to have the same 
form as that for S 2 , i.e., Eq. f )25|) . Therefore the solutions for S n (n > 3) are S n (x,t^j — ot n — 
constant. 



From (1101) . the general solution W(x,t) to FPE in the case U\ = [ix p t q , where q = —p/2, 

is 

W(x, t) cx expj - ^_ - ^ (^=Y + i (\ ai + X 2 a 2 + X 3 a 3 + ■ ■ ■ ) } 



cx exp 



In the last step, the constants a±, a 2 , are set to zero since the part exp{-^(A«i + 
A 2 «2 + A 3 a3 + ■■•)} can be absorbed into the normalization constant, which however can 
not be determined until the parameters p and /i are given. Furthermore, for W(x,t) to be 
normalizable, p can only take the values p = 0,1, 2, and p = even integer for p > 2. 

Below we shall consider two special but interesting cases of such a form of drift potentials. 



5.1 U\ = x/yji, p = 1 



case 



Let us now consider the case with \x = 1 and p = 1. The FPE is 



dW _ d^W_ \_dW_ 

~dt ~ ~dx^ + Vt~dx~ ' ( } 

Its similarity solution is 

This result is the same as that in [16]. One can see from ( )29|) that position x is shifted 
as time t elapses. 

5.2 Ui = x 2 /2t, p = 2 case 

The second example is the case with U\{x,t) = x 2 /2t, where the parameters taken are \x = 
1/2 and p — 2. From the discussions given before, we find that the probability distribution 
function 
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is the solution to the FPE 



dW 



D 



d 2 W Xx dW X 



+ -W . 



(31) 



dt 



dx 2 t dx t 



In (1301) . one can see that the perturbation parameter A controls the half width of the Gaussian 
distribution W(x,t). 

6 Summary 

When a partial differential equation possesses certain scaling behavior, the so-called similar- 
ity method is of great help in finding its solutions. One advantage of the similarity method is 
to reduce the partial differential equation to an ordinary differential equation through some 
new independent variables (called similarity variables), which are certain combinations of 
the old independent variables. The solutions so obtained, called similarity solutions of the 
differential equations, likewise possess proper scaling forms. A well-known example is the 
diffusion equation. 

In this paper, we have applied the similarity method to a class of perturbative FPE's 
with small time-dependent drift potentials studied in [16]. We have presented the main ideas 
of the similarity method. The method was then applied to find similarity solutions of the 
FPE with certain scale-invariant drift potentials. Our results show that similarity method 
can be a useful tool to solve FPE with time- dependent drift and diffusion coefficients. While 
the present work is only concerned with the perturbative FPE, it can be extended to the 
general FPE without much difficulty [19]. 
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